Rossby waves in rapidly rotating Bose-Einstein condensates 
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We predict and describe a new collective mode in rotating Bose-Einstein condensates, which is 
very similar to the Rossby waves in geophysics. In the regime of fast rotation, the Coriolis force 
dominates the dynamics and acts as a restoring force for acoustic-drift waves along the condensate. 
We derive a nonlinear equation that includes the effects of both the zero-point pressure and the 
anharmonicity of the trap. It is shown that such waves have negative phase speed, propagating in 
the opposite sense of the rotation. We discuss different equilibrium configurations and compare with 
those resulting from the Thomas-Fermi approximation. 



The rotation of Bose-Einstein condensates (BEC) has 
attracted much attention recently, both theoretically and 
experimentally [H |2] . Due to the superfluid character of 
the BEC, the effects of the rotation are quite different 
from those observed in a normal fluid and its properties 
strongly depend on the effects of the confinement. The 
pioneer experiments based in both the phase imprinting 
[3] and the rotating leaser beam [H [5] techniques inde- 
pendently confirmed the nucleation of quantized vortices, 
which is a clear manifestation of the superfluid properties 
of the condensate. Since then, much effort has been made 
to understand the dynamics of the rotating BEC [S] and, 
in particular, the mechanisms of vortex nucleation [71 [8]. 
Particular interesting features of quantized vortices in 
Bose-Einstein condensates of alkali atoms are related to 
the formation of vortex arrays, where singly quantized 
vortices typically arrange in highly regular triangular lat- 
tices, similar to the Abrikosov lattice for supercondutors. 
Such a configuration is only possible when a sufficient 
amount of angular momentum is effectively transferred 
to the system, corresponding to a situation of a rapid 
rotation. The acquired angular velocity then tends to 
enlarge the rotating cloud and the centrifugal force is re- 
sponsible for the flattening of the density proflle towards 
a two-dimensional conflguration. In the limit where the 
rotation frequency fl approaches the transverse trapping 
frequency uj± , the quadratic centrifugal and the harmonic 
trapping potentials cancel out and the system is no longer 
bounded. The possibility of reaching high angular veloc- 
ities is therefore provided by the addition of anharmonic 
terms to the trapping potential, making worthy the in- 
vestigation of new equilibrium conflgurations with differ- 
ent vortex states and new collective modes [5]. In the 
present work, we take advantage of such an interesting 
medium to predict a different hydrodynamical mode in 
rotating BECs, in complete analogy to the Rossby waves 
observed in geophysics. Rossby waves, also known as 
planetary waves, have been recognized for a long time as 
the main pattern of long period variability in the upper 
tropospheric winds [10]. These are responsible for the 
well known cyclonic and anticyclonic systems that char- 
acterize the day-to-day weather systems in mid latitudes. 



and can be observed both in the upper troposphere and 
in the oceans. The waves exist due to the variation of 
the Coriolis parameter / with latitude, which acts like a 
restoring force for an air particle that is disturbed from 
its equilibrium latitude. In a rotating BEC, the Corio- 
lis parameter is replaced by twice the angular rotation 
frequency, 2f2. We show that such waves are dispersive, 
their phase speed being always negative, which means 
that these oscillations always propagate westward with 
respect to the BEC rotation. We here derive a new equa- 
tion for the Rossby waves in a rotating condensate, which 
accounts for the vortex lattice and the anharmonicity of 
the trap. 

In the presence of a large number of vortices, Cozzini 
and Stringari showed that it is possible to average the 
velocity fleld over regions containing many vortex lines 
and assume that the vorticity is spread continuosly in the 
superfluid. This approximation is known as the diffused 
vorticity approach |llj and simply corresponds to assume 
a rigid-body rotation v = x r, where the angular ve- 
locity is $7 = flz with n = nhny/m, nu(r) is the average 
vortex density in the vicinity of r and m stands for the 
atomic mass. Therefore, the usual irrotationality condi- 
tion V X V = is no longer valid and should be replaced 
by V X v(r) — 2S},. In that case, the macroscopic dy- 
namics of the rotating BEC is provided by the rotational 
hydrodynamical equations in the rotating frame 

— + V-[n(v-f2xr)]=0 (1) 
ot 

(2) 

where v- Vv = V('y^)/2 — v x (O x v). The usual hydro- 
dynamical calculations are based on the Thomas-Fermi 
approximation, which consists in neglecting the quantum 
pressure proportional to fi^ . In this work, however, we 
include this quantum term, since we may be interested 
in Bogoliubov-like waves. This allows one to cast the ef- 
fects of the quantum depletion in the condensate, which 
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may be relevant to describe the so-called quantum tur- 
bulence [ig. Here, V(r,n) = Vt™p(r) - mfi^VS, with 
r = {x^ + y'^y^'^, stands for the effective trapping poten- 
tial, which reads 



where the material derivative can be approximated as 
d/dt « d/dt+SvQ-'V . Using the fact that In n w lnno+(f>, 
where = 6n/nao, and putting Eqs. Q, ([9| and (10) 
together, one should obtain 



(3) 



with ajio = ^/fi/ 'mLo± being the characteristic harmonic 
oscillator length and /? the dimensionless anharmonicity 
parameter. The term 2fl x v in Eq. ^ represents the 
Coriolis force, which will act as the restoring force for the 
oscillations considered here. We consider perturbations 
around the equilibirum configuration, making n = no+Sn 
and V = 5v. In that case, the system can be described 
by the following set of perturbed equations 



d_ 
di 



5n + V ■ (7io(5v) = 



(4) 
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where rioo is the peak density. The rotational velocity 
field can be split into to parts, (5v « Svq + dvp, where 



S.o - X S 



(6) 



is the zeroth-order drift velocity, resulting from taking 
d/dt = d/dt + • V in Eq. ^ and 
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(7) 



The polarization velocity (5vp is the first-order correction 
to the drift velocity ([t]) and satisfies the following equa- 
tion 



dt 



+ (5vo • V (5vo = -20 X 5wp 



(8) 



which yields 



where Sj^ — {Sx,Sy) is the transverse component of S. 
The continuity equation Q can be written in the follow- 
ing fashion 



dt 



Inn V • = 0, 



(10) 



1 - rlVl + ^+2n {V-, ^ - V^V + Inno} = 0. 

(11) 

This equation is formally similar and generalizes Char- 
ney's equation [T3], also refereed in the literature as 
Charney-Hasegawa-Mima (CHM) equation. Here, vq = 
Cs/2Q. represents the generalized Rossby radius, Cs = 
\/ gnoo/m is the sound speed, ^ = %/ ^/2mgnao is the 
healing length [Mj and = - ro^^V^(/i/2. The oper- 
ator {a, b} = r~^{dradgb — drbdga) is simply the Poisson 
bracket in polar coordinates and 9 represents the angular 
coordinate. The latter equation describes hydrodynam- 
ical drift waves in a rapidly rotating Bose-Einstein con- 
densate and includes new features relatively to the CHM 
equation, widely used in the study of the dynamics of 
waves and turbulence in plasmas and in the atmosphere. 
Namely, the terms proportional to cast the effects of 
the zero-point pressure, which are known to play no role 
in geophysics. According to typical experimental condi- 
tions, we estimate the sound speed to be ~ 1 mm/s, 
il ^ uj± sa 27r X 65 Hz [T3], a transverse harmonic oscil- 
lator radius of Uho ~ 1.7 /Ltm and a Rossby radius around 
To ^ 1.2 /im. The Rossby number, Ro, defines the ratio 
of the inertial to Coriolis forces 



Ro 



Cs 



(12) 



where L is a typical length of the system and / is the 
Coriolis parameter [101 US] • It characterizes the impor- 
tance of Coriolis accelerations arising from planetary ro- 
tation and typically ranges from 0.1 ^ 1 in the case of 
large-scale low-pressure atmospheres to 10^ in the case 
of tornados. Making L — Uho in a BEC, and / — 2n, we 
have Ro — ro/a^o ^ 0.7, which suggests that the BEC 
can be regarded as a low-pressure atmosphere. 

We now show that a rotating BEC can sustain a new 
hydrodynamic mode corresponding to a drift-acoustic 
wave. For that purpose, we keep only the linear terms 
in Eq. (Ill, and look for perturbations of the form 
(f> ^ g'(k r^:J7)^ rpj^g respective dispersion relation is then 
readily obtained and reads 



-VRkg 



l+eP/2 



H-r2/c2( 1 + ^2^/2)^ 



(13) 



where kg —\<i-eg is the polar (or zonal) component of the 



wave vector k = (/cj., ky). The term vji 



2ilrQdr In tiq 
Because the 



is the generalized Rossby (drift) velocity, 
equilibrium profile is generally very smooth, we expect vn 



3 



to be small (compared to the Bogoliubov speed Cg), which 
suggest these waves to appear as a low frequency oscilla- 
tion (compared to the both uj± and ft). The dispersion 



relation ( 13 ) is similar to the expression for barotropic 



Rossby waves in the atmosphere [TD] and to the disper- 
sion relation obtained for drift waves in a magnetized 
plasma [H]. For long wavelengths rf^k'^ ^ 1 (and con- 



sequently ^ k <^ 1), Eq. (13) reduces to the zonal 



flow dispersion relation uj w —kgvu. One of the remark- 
able features of the zonal, transverse acoustic, waves is 
that of having negative zonal phase and group veloci- 
ties, c'g^^ = c'f^ « —vr- It means that they propa- 
gate always "weastward" comparatively to the rotation 
of the condensate (which explains the negative values for 
the frequency in Eq.(13)). For short wave lengths, one 
obtains the dispersion relation for the (actual) Rossby 
waves u! « —vakg/ (rgfc^), with phase and group veloci- 
ties approximately given by 



J9) 
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In Fig. Q, we plot the dispersion relation (131 for dif- 
ferent values of the healing length ^, using vr_ — O.lcs- 
Although a single Rossby wave of arbitrary amplitude 



is a solution of the linearized wave equation Eq. (13), a 



superposition of waves, generally, is not. The nonlinear 
interaction between the waves leads to a mechanism of 
energy transference. To study the interaction properties, 
one decomposes the solution into its Fourier series, (?V = 
^jj.0kexp(ik • r) which, after plugging into Eq. (Ill, 
yields the following nonlinear equation 



d(t>k 
dt 
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where 



^ 2r25(ki+k2-k)(k2 xki).JI 
{l + ekl/2) {l + rlkl+riekt/2) 



1 



riek^/2 



(16) 



is the nonlinear coupling operator and ujk is given by 
Eq. (13). Remarkably, only the waves that satisfy the 
condition ki -|-k2 = k interact nonlinearly. The set of the 
waves satisfying this condition is known in the literature 
as the resonant triad. This resonance mechanism is able 
to transfer energy between different length scales, being 
one of the sources of classical turbulence in plasmas and 
in the atmosphere [HI [17]. Here, due to the existence of 
additional terms that properly account for the quantum 
hydrodynamical features of the system, i.e, when large 
variations of the density profile are present, we believe 



that Eq. (15) may be used to describe turbulence in 



rotating Bose-Einstein condensates, opening a stage to 
explore the similarities and differences between classical 
and quantum turbulence. 

Another interesting feature of the Rossby waves in 
Bose-Einstein condensate is the possibility of finding lo- 
calized structures, which may result, for example, from 
the saturation of the triad resonance mechanism men- 
tioned above. Such purely nonlinear solitary structures 
can be obtained from the stationary solutions of Eq. |11[ 
which readily yields 



1 



2 



= 0. 



(17) 



In the Thomas-Fermi limit, the latter expression simply 
reduces to {0, V^^^} = 0, which is satisfied for a family 
of functions — J-'{(f>), where J^{x) is an arbitrary 

function of its argument. The different choices for T will 
lead to different structures, which describe many physi- 
cal nonlinear stationary solutions. For example, for the 
choice T{<j)) oc exp{—2(f>), Stuart [18] showed that the 
so-called " cat-eye" solution describes a vortex chain in a 
magnetized plasma sheet, which has been observed ex- 
perimentally in mixing layer experiments [19 . However, 
in the present case, there are physical limitations that 
impose specific constraints to the choice of the solutions. 
In particular, as discussed in [9 , the equilibrium profile 
associated with the potential in Eq. ([3]), which is given 



by no{r) = UryoiH 



-){r^ — R_) (the peak density is 



noo — PJT-i^j_/2g), must vanish at the Thomas-Fermi radii 
defined as follows 



Rl 



"ho 



± 



2/3wj^ 



2^ 



(18) 



where ^ represents the chemical potential. For /i > 0, 
the radius i?_ is purely imaginary and the density van- 
ishes at i? = while for < both i?_ and i?+ are 
present. This reflects the transition occurring at /i = 0, 
where a hole forms in the centre of the condensate and 
the equilibrium profile assumes an annular shape. The 
simplest nonlinear structure that verifies such constraints 
is obtained for J-{(j)) = —i^4>, and the respective radial so- 
lution, for ^ = 0, yields 0(r) = AJoinr) + BYo[Kr). The 
values oi K,, A and B are such that the solution vanishes 
at the radii R±. In Fig. ([2|, we plot two possible solitary 
structures in the overcritical rotation regime 17 > wj^, 
obtained for both ^ > and fj, < 0. It is interesting to 
observe that, even for the same set of parameters, the 
resulting solitary structures may differ from the usual 
Thomas-Fermi equilibrium profiles discussed above. 

In the spirit of a mean field description, we derived an 
equation that governs the dynamics of a new drift mode 
in rotating Bose-Einstein condensates, which is in close 
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FIG. 1: (Color online) Dispersion relation of the Rossby waves 
in a rapidly rotating BEC, for vr = O.lCs. It is clearly shown 
that phase speed is generally negative. The blue full line 
corresponds to the Thomas-Fermi case, ,^ = 0. The black 
dashed and red full lines respectively correspond to ^ = 0.7ro 
and ^ = 1.3ro. 




analogy to the Rossby waves in geophysics. Our equa- 
tion, which has been estabhshed for the first time in the 
context of superfluids (for the best of our knowledge), 
casts the effect of the anharmonicity of the trap, and can 
thus be extended to the overcritical rotating regime. Af- 
ter hnearization, we derived the dispersion relation for 
the Rossby waves and showed that they propagate on 
the opposite sense of that defined by the angular rotat- 
ing frequency Vl. A particular feature of these waves is 
that they resonantly decay in form of triads, which is 
a clear manifestation of a three-wave mixing mechanism 
in BEC. A very recent and interesting work by Bludov 
et al. [22] also establishes the connection between Bose- 
Einstein condensates and geophysics, where the authors 
report on the occurrence of rogue waves, a well known 
phenomenon in deep oceans. We therefore believe that 
future numerical and experimental work would reveal in- 
teresting features concerned to the nonlinear dynamics 
in these systems, with special emphasis to the issue of 
turbulence spectrum, where we believe that the numer- 
ical integration of Eq. (15) would be relevant. A more 



detailed investigation of the temporal evolution of the 
Rossby wave turbulence, as well as the stability of the 
solitary structures certainly deserve further attention. 
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